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Abstract —  This  paper  utilizes  the  wavelet  packet  transform  to  study  the  effects  of  decreasing  the  number  of  slots  per 
pole  per  phase  on  the  harmonic  behavior  of  the  working  flux  density  in  the  air  gap.  A  surface  mounted  PM  motor  is  used  as 
an  example.  The  original  design  contains  6  PM  poles  and  36  stator  slots  with  three  phase  winding  embedded.  The  rotor  is 
redesigned  to  have  4,  8,  and  12  poles  mounted  on  the  same  rotor  iron  surface  and  the  winding  is  rearranged  to  form  4,  8, 
and  12  stator  poles  respectively  with  the  original  36  slots  stator.  The  air  gap  flux  density  waveform  along  a  round  path 
close  to  teeth  surface  is  obtained  from  the  FE  solutions.  The  harmonic  behavior  of  air  gap  flux  is  studied  with  a  wavelet 
packet  transform.  It  is  found  that  an  increase  the  number  of  poles  causes  an  increase  in  the  harmonic  content  of  the  air-gap 
flux  density  and  the  accompanying  effects  on  the  proper  operation  other  components  connected  to  external  circuits. 
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I.  Introduction 

Conventional  analytical  design  and  analysis  methods  of  AC  machines  follow  the  assumption  of  having 
sinusoidal  flux  waveform  in  the  air  gap.  An  increase  in  the  number  of  poles  can  be  utilized  for  increasing  the 
machine  power  density  or  for  satisfying  the  lower  speed  demands  at  some  application  fields.  This  means  in  turn 
reduces  the  number  of  slots  per  pole  per  phase. 

Decreasing  the  slot  number  per  rotor  pole  per  phase  may  lead  to  harmonic  behavior  increase  in  the  air  gap 
flux  waveform.  Harmonics  increase  can  be  the  source  of  a  variety  of  undesirable  effects  in  the  energy  system. 
For  example,  harmonics  can  cause  signal  interference,  over  voltages,  and  circuit  breaker  failure,  as  well  as 
equipment  overheating,  malfunction,  and  failure.  Harmonics  can  cause  excessive  heating  in  system  components, 
resulting  in  shortened  life  or  failure.  Rotor  heating  and  pulsating  output  torque  caused  by  harmonics  can  result 
in  excessive  motor  heating  and  inefficiency.  In  other  words,  the  harmonic  increase  will  affect  both  the  machine 
performance  and  the  machine  control  system  performance.  Also,  it  will  have  negative  effect  on  the  other  NEPS 
components  and  operation. 

Wavelet  transform  proves  to  have  some  advantages  when  compared  to  the  classical  FFT-based  algorithms. 
One  of  the  basic  assumptions  of  the  FFT,  the  periodicity  of  the  signal,  is  not  required  and  erroneous  results  are 
avoided.  The  wavelet  methods  solve  this  problem  by  decomposing  the  measurement  signal  into  a  family  of  time 
finite  signals  through  a  filter-bank  resulting  in  multi-scale  wavelet  decomposition.  Harmonics  of  the 
fundamental  frequency,  sub-harmonics  and  inter-harmonics  are  considered,  disregarding  the  length  of  their 
occurrence  in  time. 

With  the  help  of  FE  analysis  and  Wavelet  packet,  the  effects  of  the  decrease  the  number  of  slots  per  pole  per 
phase  on  the  harmonic  behavior  of  the  air  gap  flux  density  waveform  is  investigated. 

In  order  to  rebuild  various  pole  designs,  such  as  the  4-pole,  8-pole  and  12-pole  based  on  the  original  6-pole 
machine,  the  windings  need  to  be  rearranged  to  carry  the  three-phase  current  in  the  slots.  The  FE  computation  is 
performed  on  the  various  designs  to  obtain  the  filed  distribution  and  the  flux  density  waveform  in  the  air  gap. 
Wavelet  and  wavelet  packets  are  introduced  and  utilized  to  analyze  the  flux  density  waveforms  in  the  air  gap. 
The  effect  of  the  decreasing  the  number  of  slots  per  pole  per  phase  on  the  working  flux  is  then  discussed. 
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II.  Motor  Topologies  For  Field  Analysis 
A.  Original  design  and  redesigns 

The  structure  of  the  original  surface  mounted  PM  synchronous  motor  has  6  poles  on  the  rotor  iron  and  36 
slots  in  the  stator  iron.  The  winding  is  double  layer  and  excited  with  3  phase  supply.  Slots  per  pole  per  phase 
(spp)  are  2(36/3 *6).' The  winding  arrangement  for  complete  machine  is  shown  in  Fig.  1(a).  The  magnets  are 
radially  magnetized.  Their  operating  point  is  Br  =1.08  Tesla  and  jur  =1.07. The  angular  magnet  pole  width  is  50 
mechanical  degree.  Nonlinear  BH  curves  are  used  for  both  stator  and  rotor  iron.  In  order  to  form  different 
number  of  poles,  the  stator  winding  needs  to  be  rearranged  so  as  to  form  same  number  of  poles  on  stator  as  that 
of  rotor.  Winding  arrangement  for  different  number  of  poles  is  also  shown  in  the  Fig.l  (b),(c),(d).  Depending  on 
the  number  of  poles,  spp  will  change.  As  the  number  of  poles  goes  on  increasing,  each  phase  will  occupy  less 
number  of  slots. 


8  Pole  12  Pole 

Fig.  1  Winding  arrangement  and  pole  design  varaition 

The  redesigned  rotors  are  shown  in  Fig.l.  The  angular  magnetic  pole  width  is  kept  constant  for  all  the  pole 
variations.  The  air  gap  between  permanent  magnets  is  adjusted  accordingly  to  form  a  whole  circle.  The 
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dimensions  of  the  stator  yoke,  rotor  yoke  and  stator  teeth  in  all  the  topologies  are  kept  unchanged.  All  radial  and 
axial  dimensions  for  motor  are  same. 

Obviously,  iron  part  in  the  motor  having  higher  number  of  poles  will  be  operating  at  higher  flux  density  level 
and  goes  towards  saturation.  Saturation  is  one  of  the  important  factors  causing  harmonic  increase  in  the  air  gap 
waveforms  along  with  slotting  effects. 

B.  FEM  analysis  and  considerations 

Following  formulation  is  used  for  solving  the  transient  magnetic  solution 

3  +  Vx(i /Vx2)  =  oVF  +  Vx#c  (1) 

Hc  =  coercive  magnetic  field  of  the  magnet 
V  =  Electric  potential 


Spoles 


l2poles 


Fig.  2  Surface  flux  density  distribution  for  different  pole  number  designs 
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Nonlinear  transient  FE  analysis  was  performed  for  the  field  calculation.  Parameters  used  for  the  transient 
solution  include  time  step  of  0.3125ms  corresponding  to  3.75  degrees  of  rotor  movement,  48  steps  over  180 
degrees  of  electrical  cycle.  Moving  air  gap  feature  is  utilized  for  the  transient  operation.  The  air  gap  is  divided 
into  3  layers,  one  belonging  to  stator,  one  to  the  rotor  and  middle  one  which  actually  causes  the  rotor  to  move 
with  respect  to  stator  without  changing  the  mesh. 

Considering  motor  periodicities,  only  quarter  geometry  is  modeled  for  4,  8,  12  poles  while  for  6  poles  one 
third  geometry  is  modeled.  Results  shown  in  here  are  at  0.015sec  time  instant.  Fig.  2  shows  the  flux  density 
distribution  profiles.  Fig.  3  shows  the  air  gap  flux  density  waveform  along  the  center  locus  of  the  air  gap  at  the 
same  time  instance. 

4  poles  6  poles 


Fig.  3  Flux  density  along  the  locus  at  the  center  of  air  gap 
III.  Wavelet  Theory 

A.  Introduction 

As  in  Fourier  analysis,  the  wavelet  transform  consists  in  decomposing  a  given  function  onto  a  set  of 
“building  blocks”.  However,  as  opposed  to  the  Fourier  transformation  (FT)  in  which  the  “building  blocks”  are 
the  well-known  complex  exponentials,  the  wavelet  transform  uses  the  dilated  and  translated  version  of  a 
“mother  wavelet”  which  has  convenient  properties  according  to  time/frequency  localization.  Unlike  FT  which 
gives  a  global  presentation  of  the  signal,  Wavelet  transform  provides  a  local  representation  (in  both  time  and 
frequency)  of  a  signal;  therefore  it’s  suitable  for  analyzing  a  signal  where  Time/frequency  resolution  is  needed. 

B.  Multiresolution  analysis 

A  multiresolution  analysis  of  L2(R)  is  defined  as  a  sequence  of  closed  subspaces  Vj  of  L2(R)  ,  jeZ,  with  the 
following  properties  [1-5]: 
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•  v(x)  e  Vj  <=>  v(2x)  e  Vj+1 

•  v(x)  eF#o  v(x  + 1) e  V0 

+qo  +oo 

•  u  Vj  is  dense  in  L2  (R)  and  f|  Vj  =  {0} 

j=-°°  j=~c° 

•  A  scaling  function  ^<eFo,,  with  a  non- vanishing  integral,  exists  such  that  the  collection 
{ (p{x  <^>  l)\l  g  Z},  is  a  Riesz  basis  of  Vo 

A  signal  can  be  successively  approximated  by  wavelets  with  different  scales  (multiresolution 
decomposition).  Each  step  of  the  decomposition  of  signal  corresponds  to  a  certain  resolution.  The 
decomposition  process  can  be  iterated,  with  successive  approximations  being  decomposed  in  turn,  so  that  one 
signal  is  broken  down  into  many  lower-resolution  components.  This  is  called  the  wavelet  decomposition  tree. 
For  n-level  decomposition,  there  are  n+1  possible  ways  to  decompose  or  encode  the  signal. 

The  wavelet  packet  transform  (WPT)  method  is  a  generalization  of  wavelet  decomposition.  In  wavelet  packet 

analysis,  the  details  as  well  as  the  approximations  can  be  split.  This  yields  more  than  1  different  ways  to 
encode  the  signal.  Wavelet  packet  decomposition  is  depicted  in  Figure  4  (a) 

Let  <f)(t )  and  ^(t)  be  the  scaling  function  and  the  corresponding  mother  wavelet  function  in  the  conventional 
DWT  and  define  y/°(t)  =  (j)(t),  and  y/\t)  =  y/(t), 
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Fig.  4  (a)  Wavelet  packet  decomposition  with  successive  filtering  and  down  sampling, 
(b)  Three  level  wavelet  packet  decomposition 


For  2n  sampled  waveform  the  wavelet  transform  coefficients  (WTCs)  of  a  given  function  B(t)  at  the  J  level 
and  kth  point  can  be  evaluated  obtained  by  convolving  the  sequence  dffk)  with  a  low  pass  filter  (LPF),  and 

then  downsampling  by  a  factor  of  two.  In  the  same  manner,  coefficients  (WTCs)  of  a  given  function  B(t)  at  the 
(2i+l)th  node  can  be  evaluated  obtained  by  convolving  the  sequence  dffk)  with  a  High  pass  filter  (LPF),  and 

then  downsampling  by  a  factor  of  two.  Number  of  bands  at  jth  level  is  2j .  Mathematically,  the  (WTCs)  can  be 
written  as: 

df  ( k )  =  J  B(t)y/fk  (t)dt  =  £  h(n)d‘M  (2k  -  n) 

,  ’  "  (2) 

d]M  (k)  =  j  (t)dt  =  £  gWd'^  (2k -n) 

n 

Where 

=  f=y/li+ldfP)  =  J^g(n)V/'  Jt) 

V  2J  2 
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They  are  the  wavelet  basis  function  where  j=0, 1 , . . .  .,N  and  i  =  0, 1 , . ,  2J  1  - 1 

Since  the  wavelet  basis  are  orthogonal,  the  following  properties  should  be  fulfilled: 


jy/j,k(t)dt  =  o 

J 

J(^«)2  =  i 

J  hAWA*) = 0 


p=integer 

p=q 

p*q 

p  *  0 


(3) 


A  three  level  wavelet  packet  decomposition  tree  is  shown  in  Fig.4b.Each  tree  represents  the  signal 
decomposition  as  depicted  in  Fig.  4a. 


IV.  Wavelet  Representation  of  RMS  and  Total  Harmonic  Distortion 

In  the  wavelet  theory,  any  waveform  can  be  expressed  in  term  of  it’s  weighted  sums  of  wavelet  basis 
function  [1,2,4,6-12],  hence  B(t)  can  be  expressed  as: 

m  =  z‘  E' ‘  A  (AwZ  (t)  +  z‘  Z"  ‘  d)M  (kty™  a) = (o+ E  Z' ‘  Aiw^t)  (4) 

i- 0  k= 0  i-0  k= 0  k= 0  i=\  k= 0 


Where  d ,0  (k)  is  the  scaling  function  coefficients  . 

RMS  of  a  waveform  with  period  T  can  be  expressed  in  terms  of  wavelet  at  a  certain  level  J  as  follows: 


j  Bitfdt  =j[2y 1  a  (*)*M  (o  +  e‘ e 1  d\ ik)y/)k  it)f  dt= J[2y 1  d)  m.,  ua  dt +J  [2y 1 2y 1  <  (^;>t  iof  dq 

k= 0  i= 1  k= 0  &=0  i=l  &=0 

^  2-7-1  -1  2iV_‘/  -1 

+2j  E  E  dAw-ikWjAOAAO 

y  i=\  k= 0 

2iv-/  2J~l  -1  2JV_/  -1  f  2J  1  -1  2JV_J  -1 

=  [  E  (<,a-))2j[^(/)]2^+E  E  «w)2J[^«r^+2  E  E 

t-n  i  z— n  ;-i  t-n 


(5) 


Using  the  wavelet  orthogonality  property,  Equation  5  becomes: 

fB(t fdt  =  A'L(d‘,(k))2  (6) 

i= 0  k= 0 

The  RMS  of  the  flux  waveform  can  now  found  by  divide  by  the  wave  period  T  and  take  the  root 
square  of  equation  6. 


B„ 


2n-\ 

E  s(«)2 


n= 0 


E  E  W)2 


n= 0  k= 0 


(7) 
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B)  = 


(dj(k))2 


(8) 


Where,  Bij  is  the  RMS  value  of  frequency  band  at  node  i.  The  Total  Harmonic  band  distortion  (THBD)  is 

defined  by  the  ratio  of  the  RMS  value  of  the  harmonic  bands  at  i  >0  (i.e.  excluding  the  lowest  band)  to  the  RMS 
value  of  the  distorted  waveform.  The  THBD  for  the  flux  density  waveform  is  given  by  the  following  relation: 


THBD  =  -E  j£(B‘)2  (9) 

Brms  V  i= 1 

V.  RESULTS 

This  technique  has  been  applied  to  the  flux  density  in  the  middle  of  the  air  gap  generated  from  the  finite 
element  model.  The  files  were  then  used  as  input  to  our  Wavelet  packets  module  to  extract  the  different  harmonic 
bands  of  the  flux  density  waveform. 

A  five  level  decomposition  was  selected  to  perform  our  study.  This  will  decompose  the  signal  into  32 
frequency  band  (25 ).  The  first  eight  are  enough  to  evaluate  the  harmonic  behaviors. 

Table  1  show  the  analysis  result  of  the  flux  waveform  for  the  first  4  frequency  band  for  4,  6,  8  and  12  poles  PM- 
machines  respectively.  From  the  analysis  it’s  clear  that  there  is  increase  in  harmonics  contents  within  the  flux 
density  with  the  number  of  poles  increase.  Fig.  5  shows  the  relation  between  number  of  poles  increase  and  the 
total  band  harmonic  distortion.  There  is  a  small  increase  in  the  THBD  from  4  to  6  poles  machine.  On  the  other 
hand  bigger  increase  is  noticed  between  6  and  8  poles,  6  and  12  poles,  and  8  and  12  poles  respectively. 

From  these  results  it’s  obvious  that  there  is  an  increase  in  the  harmonics  content  with  in  the  flux  density 
waveform  with  number  of  poles  increase. 


Fig.  5  THBD  change  with  number  of  poles  change 


VI.  Conclusion 

FEM  simulation  presents  an  effective  method  for  studying  the  influence  of  number  of  slots  per  pole  per 
phase  decrease  on  the  harmonic  content  of  the  working  flux  density.  Finite  element  analysis  allows  us  to 
determine  the  flux  density  in  the  air  gap. 
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A  wavelet  packets  technique  was  used  to  evaluate  the  harmonic  behavior  for  different  number  of  poles.  The 
property  of  WPT  analysis  shows  an  ability  to  quantify  different  types  of  disturbances.  It  also  shows  high  ability 
of  wavelets  to  extract  the  different  harmonic  components  disregarding  the  length  of  their  occurrence  in  time. 
The  harmonic  contents  tend  to  increase  with  number  of  poles  increase.  The  results  of  are  very  useful  in  the 
design  and  development  of  both  the  motor  and  drive’s  diagnosis  systems. 

As  a  future  study  we  will  consider  the  effects  of  PM  Pole  structure  and  the  rotor  surface  on  the  working  flux, 
which  is  also  one  of  the  important  factors  affecting  the  harmonic  behavior  of  working  flux. 

Table  1  Total  Harmonic  distortion  for  different  number  of  poles 


THD 

% 

RMS 

Band  0 
0-75Hz 

Band  1 
75-150Hz 

Band  2 
150-225Hz 

Band  3 
225-300Hz 

Band  4 
300-375Hz 

4 

poles 

9.2 

0.5737 

0.5713 

0.0432 

0.0212 

0.0145 

0.0103 

6 

poles 

11.63 

0.5913 

0.5873 

0.0576 

0.0280 

0.0153 

0.0137 

8 

poles 

14.4 

0.5824 

0.5766 

0.0695 

0.0331 

0.0144 

0.0168 

12 

poles 

18.95 

0.5600 

0.5498 

0.0907 

0.0446 

0.0166 

0.0219 
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